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SURFACE  WAVE  STATISTICS  AND  SPECTRA 
DURING  HIGH  SEA  STATE  CONDITIONS  IN  THE  NORTH  ATLANTIC 


1.  INTRODUCTION 

The  simplest  approach  to  remotely  measuring  surface  waves  is  by 
profiling  them  from  an  aircraft  along  a  line.  This  was  discussed  by 
Longuet-Higgins  [1956]  but  few  results  using  the  technique  were  published 
until  Barnett  and  Wilkerson  [1967].  Since  then,  several  studies  have 
been  conducted  mostly  on  wind-generated  wave  fields.  Most  notable  works 
are  by  Schule  et  al.  [1971]  and  several  by  Ross  including  a  recent  study 
on  Lake  Michigan  [Liu  and  Ross,  1980]. 

The  use  of  airborne  platforms  offers  several  advantages  and  in  the 
case  of  the  profilometer,  some  disadvantages.  Large  areas  can  be  surveyed 
in  a  short  amount  of  time  allowing  flexibility  in  sampling  schemes 
especially  in  complex  or  Inhomogeneous  wave  fields.  Perhaps  the  most 
important  advantages  of  profilometry  are  the  ease  of  deployment,  the 
reliability  of  hardware,  and  the  simplicity  of  the  analysis  as  compared 
to  most  other  devices  that  can  be  used  in  deep  water.  Examples  of 
situations  where  the  mobility  of  airborne  profilometry  is  of  definite 
advantage  is  for  sea  truthing  other  remote  sensors,  especially  spaceborne, 
and  in  wind  wave,  wave-current  and  perhaps  wave-bottom  interaction 
studies.  The  disadvantages  are  that  the  seas  rarely  consist  of  purely 
swell  or  purely  wind  driven  seas  and,  therefore,  any  one  flight  track 
will  result  in  aliasing  of  at  least  one  major  component.  This  can  be 
overcome  somewhat  by  flying  several  tracks  in  a  pattern.  Also,  in  cases 
where  sharp  changes  in  sea  state  occur,  the  spectrum  may  not  show  the 
true  surface  change  since  a  finite  distance  is  required  to  generate  an  average 
spectrum  with  a  reasonable  number  of  degrees  of  freedom  and  therefore, 
smoothing  results. 

The  data  set  used  in  this  report  was  collected  for  the  purpose 
of  validating  the  GE0S-3  radar  altimeter  significant  wave 
height  (H.,-)  estimates.  Initial  results  from  that  experiment  in  1976 
have  been  published  by  McClain  et  al.  [1979].  The  data  set  consists  of 
six  missions  flown  during  high  sea  states  and  is  of  sufficient  quality 
to  allow  the  testing  of  more  recent  theoretical  and  empirical  results 
which  imply  that  additional  information  may  be  extracted  from  altimeter 
data. 


According  to  the  theoretical  development  by  Longuet-Higgins  [1963], 
the  third  cumulant  and  therefore  the  skewness  of  the  distribution  of 
surface  elevation  can  be  evaluated  in  terms  of  the  frequency  spectrum. 
This  formalism  has  been  applied  by  Huang  and  Long  [1980]  who  incorporated 
the  analytical  form  of  the  Phillips  spectrum.  The  result  was  a  linear 
relationship  between  skewness  and  significant  slope  (rms  wave  height/ 
dominant  wavelenghc).  The  relationship  was  corroborated  by  their 
wind-wave  tank  data.  Longuet-Higgins  on  the  other  hand  applied  his 
analysis  to  field  data  by  Kinsman  [1960]  and  did  not  incorporate  any 
theoretical  models  of  the  wave  spectrum  into  his  analysis.  It  is 
important  to  note  that  measurements  of  skewness  in  the  literature  are 
rare. 
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Over  Che  past  Chircy  years,  a  great  amount  of  effort  has  been  directed 
towards  finding  an  analytical  expression  for  the  wave  spectrum.  The 
most  recently  proposed  model  is  the  Wallops  Spectral  Model  (WSM)  by 
Huang  et  al.  [1981].  This  model  requires  two  inputs,  the  rms  wave 
height  and  the  dominant  wavelength.  Like  previous  models,  it  assumes 
that  the  wave  field  has  only  one  dominant  component  and  is  not  a  superposition 
of  swell  and  wind  waves  with  comparable  energies. 

Besides  being  mathematically  simple  in  form  and  requiring  only  two 
inputs,  the  model  has  an  interesting  connection  to  satellite  radar 
altimetry  which  couples  with  the  skewness  relationship.  Walsh  [1979] 
has  shown  that  it  may  be  possible  to  measure  skewness  with  an  altimeter. 

His  results  were  derived  from  GEOS-III  data  and  although  his  results 
will  be  critically  discussed  later,  the  possibility  has  been  demonstrated. 

For  the  purposes  of  this  paper,  it  is  assumed  that  accurate  measurement 
of  skewness  from  an  altimeter  is  possible  since  an  altimeter  such  as 
the  one  on  SEASAT  provides  much  more  accurate  and  detailed  information 
on  the  return  waveform.  Also,  the  rms  wave  height  has  been  proven  to 
be  a  reliable  measurement  from  an  altimeter  [Parsons,  1979].  Thus,  by 
applying  the  skewness  relationship,  the  dominant  wavelength  can  be 
estimated  and  together  with  the  rms  waveheight,  the  ocean  wave  spectrum 
can  be  calculated  using  the  WSM. 

It  is  the  Intention  of  this  paper  to  present  results,  derived  from 
airborne  profilometer  measurements  of  surface  waves,  and  to  compare  them 
with  those  predicted  by  Huang  and  Long  [1980]  and  Huang  et  al.  [1981]. 

In  particular,  we  want  to  test  the  relationship  between  skewness  and 
significant  slope  for  a  variety  of  combinations  of  sea  and  swell. 

Secondly,  we  want  to  investigate  the  possibility  of  using  skewness 
measurements  for  calculating  parameters  in  the  Wallops  model.  Finally, 
we  want  to  check  the  spectral  model  when  the  parameters  are  calculated 
from  field  data  with  swell  effects  removed. 

2.  REVIEW  OF  THEORY 

Huang  and  Long  [1980]  have  shown  that  the  available  skewness  data 
compares  very  well  to  the  theoretical  relationship, 

K3  -  8tt£  (1) 

where  K  is  the  skewness  of  the  surface  wave  elevation  distribution 
relative  to  the  mean  sea  level  and  £  is  the  significant  wave  slope  which 
is  defined  as 


C  -  <S2)1/2A0  (2) 

2 

c,  is  defined  as  the  mean  square  surface  elevation  relative  to  the  mean 
sea  level  and  X  is  defined  as  the  dominant  wavelength  at  the  peak  of 
the  ocean  wave  spectrum.  From  (1)  the  significant  wave  slope  £  can  be 
evaluated  if  the  value  of  skewness  is  known  from  the  measurement. 
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Furthermore,  from  the  mean  square  surface  wave  amplitude  and  (2),  the 
dominant  wavelength  X  will  be  calculated.  Knowledge  of  the  value  of  £ 
and  \Q  are  essential  ?o  construct  the  WSM.  Considering,  a  single  train 
of  waves  the  slope  m  of  the  line  connecting  the  principle  and  the  first 
harmonic  waves  in  the  Stoke' s  expansion  should  be 

m  -  log(2ako)2  (3) 

log  (1/2) 

Extending  this  result  to  a  broader  band  wave  spectrum  which  is  a  better 
representation  of  the  random  wave  field,  (3)  will  have  the  form  of 

•-Kr§2i| 

since,  for  a  random  wave  field. 


C2  -  a2/ 2 


(5) 


and 

ko  *  f 
o 

Once  the  significant  wave  slope  5  is  evaluated,  from  (4),  m  is  known. 

The  Wallops  wave  spectral  model  has  the  form  of 

2 

4>(n)  *  (nQ/n)m  exp  [-<5(no/n)n]  (7) 

n  ^  u 

o 


where  n  is  the  wave  frequency  at  the  spectral  peak,  g  is  the  gravitational 
acceleration,  and  S,  6  and  n  are  coefficient  functions.  Since 
3<Kn)/3n  ■  0  at  n  ■  n  ,  it  is  concluded  that  m  ■  fin.  If  0  is  chosen  to 
be  4  as  did  Pierson  and  Moskowitz  [1964]  and  Hasselmann  et  al.  [1976], 

6  -  m/4  (8) 

Furthermore,  with  the  definition  of  the  spectral  function  as 


?2  ■  /q  <Kn)dn 

2 

and  incorporation  of  the  linear  dispersion  relationship  of  nQ 


8V 


g  -  (2TTC)2m(a~1)M 
r^m-l^  .(a-57/4 


(9) 
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Where  T  is  the  Gamma  function.  Therefore,  from  (4)  to  (9),  the  Wallops 
wave  spectral  model  is 


*(n)  *  m  '  (f?«)  exP  I'  f(tl0/n4]  (10) 

n  n 

o 

Following  the  theoretical  work  laid  out  by  Longuet-Higgins  [1963], 
the  theoretical  bounds  of  the  skewness  are 

5.28iTC(m-l)/(m-2)  <  K3  <  12.12Tr£(m-l)/(»-2)  (11) 

Equations  (1)  to  (11)  form  the  basic  equations  for  the  comparisons  between 
the  theoretical  predictions  and  the  measurements  made  during  the  experiment. 

3.  EXPERIMENT 

An  airborne  experiment  was  conducted  off  the  coast  of  Newfoundland, 

Canada,  in  February  1976  (see  Figure  1).  The  original  intent  of  the 
experiment  was  to  collect  aircraft  measurements  of  Significant  Wave 
Height  (SWH)  to  be  used  to  test  and  validate  the  performance  of  the 
GEOS- 3  satellite  altimeter  algorithms.  The  experiment  was  conducted  by 
the  Space  Sensing  Applications  Branch  at  the  Naval  Research  Laboratory  (NRL) . 

One  of  the  instruments  aboard  the  NASA  Wallops  Flight  Center  (NASA/WRC) 

C-54  aircraft  was  a  Geodolite  continuous  wave  laser.  The  significant  wave 
height  was  reported  to  be  4.0  -  8.0  meters  by  Parsons  [1979]  and  McClain 
et  al.  [1979].  For  the  present  study,  more  refined  analysis  of  the  laser 
data  resulted  in  improved  values  of  significant  wave  height.  Figure  2  is 
taken  from  McClain  et  al.  [1979]  but  with  the  new  laser  values.  The  sea  states 
during  the  six  flights  were  generally  high  and  were  composed  of  various 
combinations  at  sea  and  swell.  The  wave  profiles  thus  constitute  a 
unique  data  set  and  is  most  appropriate  for  testing  the  general  applicability 
of  the  relationships  discussed  in  section  2. 

4.  DATA  ANALYSIS 

In  order  to  measure  statistical  moments  of  higher  order  than  the  second, 
an  improved  approach  implemented  for  removing  aircraft  motion  from 
the  profilometer  record  over  that  previously  used  [McClain  et  al.,  1980] 
is  required.  The  earlier  approach  was  to  incorporate  verticle  aircraft 
acceleration  into  the  correction  scheme  by  (a)  filtering  both  laser  and 
accelerometer  records  using  the  same  filter  (b)  Fourier  transforming 
them  separately,  (c)  doubly  Integrating  the  acceleration  power  spectrum  by 
dividing  each  component  by  a  factor  of  (2irf )  where  f  is  the  frequency; 
and  (d)  subtracting  that  spectrum  from  the  profilometer  spectrum.  To  generate 
a  corrected  time  series  from  which  higher  moments  can  be  easily  calculated, 
the  double  integration  was  performed  in  the  time  domain.  Parabolic  trends  result 
and  various  techniques  have  been  devised  to  remove  them.  Several  schemes 
Including  a  least  squares  curve  fitting  to  a  parabola  were  tested. 
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but  It  was  found  that  by  first  subtracting  the  mean  acceleration  from  the  data 
segment,  then  double  integrating  and  finally  applying  a  numerical  filter, 
the  most  satisfactory  result  was  produced.  The  procedures  were  tested 
using  artifically  generated  laser  and  accelerometer  time  series. 

Accelerometer  data  is  used  in  both  frequency  and  rime  domain  processing 
to  remove  residual  vertical  aircraft  motion  that  the  filter  passes. 

The  filter  is  designed  with  a  cut-off  which  allows  reasonable  wavelengths 
of  encounter  to  pass.  The  final  wave  record  is  obtained  by  subtracting 
the  filtered  profllometer  series  from  the  final  aircraft  motion  record 
since  the  laser  measures  the  inverted  surface.  Also,  it  appears  that 
the  laser  phase  shift  correction  routine  needed  for  this  particular 
instrument  can  induce  spurious  long  wavelength  components  making  comparison 
of  long  data  records  of  laser  and  aircraft  motion  difficult.  A  pulsed 
laser  system  such  as  is  presently  used  by  NRL  measures  absolute  range 
and  does  not  encounter  this  problem.  Filtering  removes  this  error 
source. 

The  filter  used  was  a  nonrecursive  Martin  filter  [Martin,  1959]  and 
its  design  for  this  application  is  discussed  in  McClain  and  Walden 
[1979].  It  has  the  necessary  property  of  zero  phase  shift  since  the 
weighting  function  is  symmetric.  Figure  3  provides  the  response  function 
for  the  filter  applied  in  all  but  one  instance. 

Before  Fourier  analysis,  each  data  segment  to  be  transformed  was 
smoothed  using  a  4-term  Blackman-Harris  window  as  described  in  Harris 
[1978].  This  window  is  characterized  by  one  main  lobe  with  a  6  dB  band 
width  of  2.7  bins  and  no  side  lobes.  All  final  spectra  were  derived 
by  averaging  five  spectra  each  being  the  Fast  Fourier  Transform  (FFT) 
of  2048  points.  The  digitization  rate  was  45  Hertz.  This  corresponds 
to  a  distance,  D,  given  by 


D 


(2048) 

45 


(12) 


where  U  is  the  ground  velocity  of  the  C-54  aircraft  in  meters  per  second. 
The  slowest  ground  speed  occurred  for  data  set  2/23b,  i.e.,  data  set  b 
on  Feb.  23,  and  was  50  m/s.  Therefore,  D  ■  2.28  km  or  about  nine  dominant 
wavelengths  on  that  day  and  the  sample  spacing  was  about  1  meter.  It 
was  felt  that  a  transform  of  this  size  was  appropriate  in  order  to 
provide  sufficient  spectral  resolution  for  our  application,  i.e.,  the 
determination  of  dominant  wavelength. 

The  algorithm  for  true  frequency,  f is  based  on  the  Doppler  shift 
expression, 

2irf  -  2?rf  +  k  •  U  (13) 

c  a  c 

where  f  is  the  apparent  frequency,  k£  is  the  wave  number  vector  associated 
with  f  ,  and  U  is  the  ground  velocity  of  the  aircraft.  For  sea  states 
the  size  investigated  and  upwind  ground  speeds  being  rather  slow,  the 
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inclusion  of  f  in  (13)  is  necessary.  For  instance,  Che  calculated  f£ 
for  the  dominant  wave  component  on  2/23b  corresponds  to  a  phase  speed  of 
approximately  20  m/s  or  401  of  the  ground  speed  and  negates  any  assumption 
that  the  surface  is  relatively  stationary. 

The  facts  that  ground  speed  varies  from  one  flight  line  to  the 
next,  that  the  transforms  were  uniform  in  size,  and  that  (13)  is  nonlinear 
in  f  do  lead  to  variable  resolution  in  the  frequency  domain.  The 
spectra  from  2/23a  and  2/23b  exemplify  this  since  the  ground  speed  of 
2/23a  was  120  m/s. 

For  the  flights  on  2/20,  2/24,  2/28,  and  3/3,  significant  amounts 
of  swell  were  present.  The  algorithm  assumes  only  one  angle  of  encounter 
and  since  all  flights  were  up  or  down  wind,  only  the  wind  wave  components 
were  accurately  mapped.  Therefore,  the  points  labeled  swell  are  not 
accurately  located,  but  no  method  to  correct  them  is  possible  with  a 
flight  pattern  with  one  leg.  This  problem  is  particularly  bothersome  on 
2/24.  The  effect  of  sampling  a  directionally  dispersed  wind  wave  field 
with  a  line  has  been  analytically  studied  in  Hammond  and  McClain  [1980] 
and  is  shown  not  to  be  a  serious  problem.  The  moments  were  generated 
from  ungrouped  data,  i.e.,  the  probability  density  functions  were  not 
determined.  These  values  were  used  in  determining  the  cumulants,  K  . 

The  equations  for  the  Kn's  are  the  same  as  those  in  Huang  and  Long  ?1980]. 

5.  RESULTS 

The  ocean  wave  spectra  (Figures  4  through  9)  can  be  typically 
classified  into  three  types  of  sea  states.  The  first  type,  represented 
by  Figure  4  through  6  show  the  sea  state  is  not  only  wind-generated  but 
also  contains  wave  energy  from  the  swell.  The  dominant  wind  wave  peak 
has  been  identified  as  W,  and  the  dominant  swell  peak  by  S.  Generally 
speaking,  the  total  wave  energy  contained  in  the  various  swell  components 
is  not  large  in  comparison  with  the  total  wave  energy  of  the  wind  wave 
components.  The  significant  wave  height  for  2/23 a  was  7.89 
meters,  and  wave  spectrum  at  high  wave  frequencies  approached  those 
predicted  by  Phillips  [1958]  assuming  the  equilibrium  coefficient 
3  to  be  0.81vl 0  .  The  predicted  peak  value  uging  Phillips  model  and 

the  observed  value  for  2/23a  are  152  and  153  m  s  respectively.  The 
second  type  of  ocean  wave  spectra  is  represented  by  Figure  7  for  the 
date  of  2/23b.  The  significant  wave  height  was  8.7  meters  and  the  sea 
state  was  purely  wind-generated.  The  third  type,  represented  by  Figures 
8  and  9,  show  that  the  sea  is  essentially  swell. 

The  wave  statistics  and  the  hindcasts  for  all  of  these  dates  are 
presented  in  Table  1.  The  dominant  wavelengths,  X  ,  were  obtained 
separately  for  the  dominant  swell  and  wind  waves.  In  case  there  was 
inadequate  information  available,  the  space  is  left  blank.  A  long  bar 
indicates  that  no  value  was  determined.  The  values  of  H.  ,y  K y  the  fourth 
order  cumulant  K^,  the  fifth  order  cumulant  K^,  and  hindcasts  on  wind  vectors 
have  been  evaluated  or  observed  for  the  combined  sea;  therefore,  their 
values  are  not  particularly  associated  with  either  swell  or  wind  waves. 
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The  significant  slopes  E  as  defined  by  (2)  were  calculated  separately  for 
wind  waves  and  swell  by  using  the  combined  rms  wave  amplitude  shown  in  Column  3, 
and  the  dominant  wavelengths  for  wind  waves  and  swell  as  shown  in  Column  4. 

The  values  of  X  /H.  ..  shown  in  Column  9  are  one  quarter  of  the  inverse 
values  of  Various  hindcast  values  of  the  wave  field  and  winds  are  in 

Columns  10  and  11,  respectively,  and  were  taken  from  Mennella  et  al. 

(unpublished  manuscript,  1976).  In  column  12,  the  equilibrium  range 
slopes  are  estimated  from  Figures  4  through  9. 

Figure  10  shows  the  plot  of  K^  with  respect  to  £.  The  data 
provided  in  Huang  and  Long  [1980]  are  shown  as  small  solid  circles 
(wave  tank)  and  A  (field  data  from  Kinsman) .  The  solid  curve  represents 
the  empirical  relationship,  (1).  On  the  dates  2/28  and  3/3  two  data  points 
are  shown.  The  data  on  the  left  indicates  swell  peaks  and  the  data 
on  the  right  indicates  wind  wave  peaks.  They  are  connected  by  a  horizontal 
line  whose  length  is  an  indication  of  the  range  of  variation  of  £.  As  one 
can  observe  from  Figure  10  the  data  points  of  the  present  study  are  coincident 
with  those  of  Kinsman  but  have  less  scatter.  The  wavelength  used  for  2/24 
is  the  hindcast  value,  since  the  flight  track  was  parallel  to  the  wind  and 
apparently  at  a  large  angle  to  the  swell.  A  special  filter  had  to  be  applied 
in  order  to  pass  the  aliases  swell  component. 

At  this  point  a  comment  regarding  the  results  of  Walsh  [1979]  is 
appropriate.  The  GEOS-3  data  represented  in  that  paper  corresponds  to 
our  underflight  data  sets  2/20  and  2/23a.  The  skewness  values  determined 
by  Walsh  from  the  GEOS-3  altimeter  were  0.21  and  0.54,  respectively, 
as  compared  to  our  values  of  0.11  and  0.20  derived  from  the  profilometer . 

In  order  to  evaluate  £,  the  hindcast  wavelengths  used  by  Walsh  were  150  m 
for  2/20  and  96  m  for  2/23.  Although  we  cannot  determine  the  swell 
wavelength  on  2/20  precisely,  the  wind  wave  peak  was  96  m  and  on  2/23a 
the  peak  was  242  m.  The  2/20  data  was  swell  dominated.  Therefore,  it 
appears  that  the  apparent  agreement  between  Walsh's  results  and  (1) 
occurred  because  the  consistently  over  estimated  K^  was  balanced  by 
under  estimation  of  X  .  We  find  that  the  hindcast  estimates  of  SWH 
(Figure  2)  and  of  wavelength  of  period  (Table  1)  tend  to  be  low  at 
least  for  these  data  sets. 

Figure  11  shows  the  plot  of  the  fourth  order  cumulent  K^  with 
respect  to  £.  Again  the  data  points  of  the  present  study  stay  within 
the  same  range  of  E  values  as  those  of  the  other  in  situ  field 
study  but  with  less  scatter.  Figure  12  shows  the  plot  of  the  fifth  order 
cumulent  with  respect  to  the  significant  slope.  The  solid  curve  is  the 
empirical  relationship. 


k5  =■  -110  E  (14) 

by  Huang  and  Long  [1980].  It  should  be  noted,  however,  there  is  no 

in  situ  field  data  available.  The  values  are  well-behaved  and  comparatively 

uniform  in  value. 

In  order  to  illustrate  more  effectively  the  skewness  dependence,  a 
linear  plot  with  respect  to  X  /H.  (Column  9  in  Table  1)  is  shown  in 
Figure  13.  The  dashed  curves  are' the  theoretical  upper  and  lower  bounds. 
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as  shown  by  (11).  The  A  points  are  defined  later.  Because  various 
proportions  of  the  total  wave  energy  are  contributed  by  swell,  the  data 
points  for  the  date  of  3/3  are  outside  the  bounds.  The  data  points  of 
2/24  are  not  included  here.  Otherwise,  the  data  points  for  the  remaining 
dates  are  reasonably  predicted  by  (1),  and  the  theoretical  bounds. 

Since  the  significant  wave  height  used  includes  both  the  swell  and  wind 
wave  contributions,  it  is  indeed  arguable  whether  or  not  this  is  a  fair 
comparison.  Longuet-Higgins'  [1963]  original  derivation  should  be  applicable 
for  the  combined  sea  state,  but  the  theoretical  curves  were  derived  by 
using  a  continuous  wave  spectrum  extending  from  the  dominant  wave  through 
the  high  frequency  wave  domain. 

If  skewness  and  rms  waveheight  are  the  measured  quantities,  all  of 
the  terms  in  the  Wallops  spectral  model  can  be  calculated  and  compared 
with  the  values  in  Table  1.  This  procedure  simulates  an  altimeter- 
derived  product.  Table  2  shows  the  parameters  derived  from  (1),  (2), 

(4),  and  (9)  using  the  values  of  skewness  and  of  Table  1.  The 

spectra  shown  in  Figure  15  are  two  modeled  spectra  and  a  profilometer 
spectrum.  The  model  spectrum  derived  from  direct  evaluation  of  dominant 
wavelength  and  rms  wave  height  shows  excellent  agreement.  The  second 
modeled  spectrum  is  derived  using  the  measured  and  (1).  The  fit  is 
not  bad  although  the  dominant  wavelength  is  now  only  198  m. 

Of  particular  interest  to  the  current  study  is  the  consideration  of 

cases  where  the  effects  due  to  swell  on  H,  ,  5,  m,  3,  A  ,  n  ,  and  f 

l/3oo  o 

are  removed.  From  the  ocean  spectra  it  is  easy  to  remove  the  swell 
contributions  to  the  root  mean  squared  surface  wave  amplitudes.  The 
arrows  in  Figures  4  to  9  indicate  the  separation  of  sea  and  swell.  The 
separation  point  was  determined  by  the  slope  of  the  equilibrium  range 
and  by  using  the  hindcast  predictions  of  seas  and  swell.  The  result  of 
this  separation  is  given  as  E  /E,  the  ratio  of  swell  energy  to  total 
energy  and  Is  listed  in  Table  1.  The  relationship  of  significant  wind 
wave  height  to  this  quantity  is 


H 


w 


1/3 


(1  E  }  Hl/3 


(15) 


Using  the  observed  wind  wave  peaks,  n  is  known  and  Table  3  can  be 
generated.  The  superscript  W  indicates  Wind  Waves. 

The  skewness  contribution  due  to  swell  is  difficult  to  evaluate  and 
cannot  be  easily  separated  in  these  data  since  the  swell  is  mismapped. 
Such  a  separation  would  assume  that  the  swell  does  not  interact  with  the 
wind-wave  field.  Perhaps  a  calculation  of  the  bispectrum  (Hasselmann, 
et  al,  1973)  would  be  useful.  In  order  to  evaluate  the  effect  of  swell 
in  a  quantitative  manner,  H.  . _ / X  is  calculated  from  Table  3  and  plotted 
as  A  on  Figure  13.  Using  (I;  or°the  solid  curve  in  Figure  13,  the 
predicted  values  of  skewness  K’^  is  determined.  Figure  15  is  a  plot  of 
the  difference  between  the  measured  and  the  predicted  skewness,  K-  -  K',, 
as  a  function  of  E  /E.  A  trend  seems  to  exist.  Figure  16  presents  the 
standard  deviation  of  the  estimates  for  each  data  set  as  a  function 
of  Eg/E.  As  expected,  the  variabilities  in  are  relatively  high 


compared  to  the  mean  values  but  the  standard  deviations  are  similar  to 
those  of  the  tank  data  in  Figure  10.  No  apparent  trend  is  seen  in 
Figure  16. 

Finally,  using  the  parameters  in  Table  3,  a  comparison  of  the  wind-wave 
portion  of  3/3  and  the  model  is  shown  in  Figure  17.  The  fit  is  quite 
good  and  illustrates  the  point  that  subfields  can  be  modeled  independently 
and  the  wave-wave  interaction  between  swell  and  the  wind-wave  field  are 
weak. 


6.  CONCLUSIONS 

Regarding  skewness,  the  data,  though  limited,  tends  to  support  the 
skewness  relation,  (1),  even  with  moderate  amounts  of  swell  present. 
However,  as  the  swell  contribution  to  total  energy  increases,  larger 
deviations  result.  This  is  expected  since  swell  should  be  more  Gaussian. 
Computation  of  and  help  demonstrate  the  reliability  of  these 
measurements  by  comparing  favorably  with  other  data  and  exhibiting 
values  appropriate  to  nearly  Gaussian  processes.  More  data  from  swell- 
dominated  seas  needs  to  be  analyzed  and  large  quantities  probably 
exist  but  are  left  unstudied  since  most  interest  has  been  in  wind  wave 
generation. 

Comparisons  of  the  WSM  and  measured  spectra  show  excellent  agreement 
and  indicate  that  sub-fields  can  be  modeled  independently.  Also,  the 
feasibility  of  using  the  skewness  and  rms  height  measurements  from  an 
altimeter  to  infer  wave  spectra  using  the  WSM  is  promising,  especially 
in  wind-wave  dominated  seas.  These  seas  are  of  greater  practical 
interest  to  mariners  and  for  wave  forecasting.  It  remains  to  be  proven, 
however,  that  the  altimeter  can  accurately  measure  skewness. 
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Table  1  -  Wave  Statistics  and  Hindcasts  from  Prof llometer .  T,  Wave  Period  in 
Seconds;  E  ,  Wave  Energy  Due  to  Svell;  E,  Total  Wave  Energy 


Table  2  -  Parameters  for  Probability  Distributions  and  Spectra 


DATE 

SEA  STATE 

Eq.  (1) 

K 

Eq.  (4) 
m 

Eq.  (9) 

B 

x  10* 

Eq.  (2) 

\ 

o 

n 

o 

f 

o 

2/20 

SWELL 

AND 

WIND  WAVE 

0.00437 

11.37 

3.18 

271 

0.4767 

0.0759 

2/23n 

SWE1.L 

AND 

WTND  WAVE 

0.00788 

9.67 

6.14 

250 

0.4963 

0.0790 

2/23b 

SWELL 

AND 

WTND  WAVE 

0.01095 

8.72 

8.76 

198 

0.5577 

0.0888 

2/24 

NO  El 

[TRIES 

2/2B 

SWELL 

AND 

WIND  WAVE 

0.00532 

10.81 

3.98 

210 

0.5415 

0.0862 

3/3 

SWELL 

AND 

WIND  WAVE 

0.00318 

12.3 

2.23 

433 

0.4490 

0.0714 

Table  3  -  Revised  Parameters  for  Wind  Wave  Spectra  by 
Removing  the  Effects  of  Swell  on  Sea  State 
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NEWFOUNDLAND  SEA  TRUTH  MISSION 
February -March  1976 


SEA  TRUTH  Hl/3(Meters) 

Fig.  2  -  Comparison  of  GEOS-3  and  Sea  Truth  Significant  Wave  Heights 


f 


r 


I  I  I  l  i  i  i 

5.6  2.8  1.9  1.4  rQ  1.1  0.9 

NORMALIZED  FREQUENCY,  r  x  10s 


Fig.  3  -  Response  Function  for  the  Filter 
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UOQ  (SPECTRA) 
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Fig.  4  -  Ocean  Wave  Spectrum,  Feb.  23a,  1976.  W,  Wind  Wave  Peak;  S 
Swell  Peak;  +  ,  Separation  of  Swell  and  Wind  Wave 
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LOG  (SPECTRA) 


The  Plot  of  Skewness  K.  with  Significant  Slope  e.  Data  of 
the  Wave  Tank  Study  (Huang  and  Long,  1980) ,  with  One-Sigma 
Errors  Bars;  A,  Results  of  In-Situ  Field  Study  (Kinsman,  1960) 


The  Plot  of  Fourth-Order  Cumulant  with  Significant  Slope  £ 
For  legend,  see  Fig.  10 
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Fig.  12  -  The  Plot  of  Fifth-Order  Cumulant  with  Significant  Slope 

For  legend,  see  Figure  10. 
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K,  and  X  /H.  with  Swell,  Wind-Wave,  and  Revised 

Wind-Wave  P efcKs  for  the  Dates  of  2/20,  2/23a,  2/23b,  2/28,  and  3/3 


Power  Spectral  Density 


